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Abstract 

This chapter reviews the theoretical description of interacting fermions in one 
dimension. The Luttinger liquid concept is elucidated using the Tomonaga- 
Luttinger model as well as integrable lattice models. Weakly coupled chains 
and the attempts to experimentally verify the theoretical predictions are dis- 
cussed. 
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1. Introduction 

In this chapter we attempt a simple selfcontained introduction to the main 
ideas and important computational tools for the description of interacting fermi- 
ons in one spatial dimension. The reader is expected to have some knowledge 
of the method of second quantization. As in section 3 we describe a con- 
structive approach to the important concept of bosonization, no quantum field- 
theoretical background is required. After mainly focusing on the Tomonaga- 
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Luttinger model in sections 2 and 3 we present results for integrable lattice 
models in section 4. In order to make contact to more realistic systems the 
coupling of strictly Id systems as well as to the surrounding is addressed in sec- 
tion 5. The attempts to experimentally verify typical Luttinger liquid features 
like anomalous power laws in various correlation functions are only shortly 
discussed as this is treated in other chapters of this book. 

2. Luttinger liquids - a short history of the ideas 

As an introduction the basic steps towards the general concept of Luttinger 
liquids are presented in historical order. In this exposition the ideas are dis- 
cussed without presenting all technical details. This is done in section 3 by 
disregarding the historical aspects aiming at a simple presentation of the im- 
portant practical concepts like the " bosonization of field operators". 



In a paper on incoherent x-ray diffraction Bloch [1] realized and used the 
fact that one-dimensional (d = 1) noninteracting fermions have the same type 
of low energy excitations as a harmonic chain. The following discussion of 
this connection is very different from Bloch's own presentation. 
The low energy excitations determine e.g. the low temperature specific heat. 
Debye's famous T 3 -law for the lattice contribution of three dimensional solids 
reads in d = 1 



where c s is the sound velocity. At low temperatures the electronic contribution 
to the specific heat in the "Fermi gas" approximation of Pauli is also linear in 
T and involves the density of states of the noninteraction electrons at the Fermi 
energy. This yields for spinless fermions in d = 1 



where vp is the Fermi velocity. With the replacement c s <-> vp the results are 
identical. This suggests that apart from a scale factor the (low energy) exci- 
tation energies and the degeneracies in the two types of systems are identical. 
For the harmonic chain the excited states are classified by the numbers n k of 
phonons in the modes uj k whith n k taking integer values from zero to infinity. 
The excitation energy is given by E({n km }) - E = J2k m nuJ k m n km - For 
small wave numbers k m the dispersion is linear uj km ~ c s \k m \. Therefore the 
excitations energies are multiples of hc s (2ir/L) for periodic boundary con- 
ditions and multiples of A# = fic s ir/L for fixed boundary conditions. The 
calculation of the partition function is standard textbook material. This is also 
true for noninteracting electrons but there the calculation involves fermionic 
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Bloch's method of "sound waves" (1934) 




(l.l) 




(1.2) 
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occupation numbers which take values zero and one. The two textbook cal- 
culations yield Eqs. (1.1) and (1.2), but through the "clever" use of the grand 
canonical ensemble in order to simplify the fermionic calculation the identity 
(apart from c s <-> vp) remains mysterious. A deeper understanding involves 
two steps: 

1) Linearization of the kinetic energy = Ti 2 k 2 / (2m) of the free fermions 
around the Fermi point kp for fixed boundary conditions or both Fermi points 
±kp for periodic boundary conditions. As the argument is simplest for fixed 
boundary conditions [2] which lead to k m = mir/L we discuss this case for the 
moment. Then the energies Ek n — £f are integer multiples of A^ = hvpir/L 
where vp is the Fermi velocity. 

2) Classification of any state of the Fermi system by the number rij of up- 
ward shifts by j units of Ap with respect to the ground state. As the fermions 
are indistinguishable the construction of the {rij} shown in Fig. 1.1, where 
the highest occupied level in the excited state is connected with the highest 
occupied level in the ground state and so forth for the second, third . . . highest 
levels, completely specifies the excited state. As the rij can run from zero to 
infinity like bosonic quantum numbers and the excitation energy is given by 
^Zj{jAp)nj the canonical partition function for the noninteracting fermions 
has the same form as the canonical partition function for the harmonic chain 
apart from A^ <-> A# if one assumes the Fermi sea to be infinitely deep [3]. 
As we have linearized u)k for small k as well as Ek around kp this equivalence 
only holds for the low temperature specific heats (&bT <C hco max , fce^ < 

Ep). 

If we denote the creation (annihilation) operator of a fermion with k n = 
nir/Lby c]j(c n ) and assume a strictly linear dispersion = fivpk n for all 
k n > a more technical formulation of the discussed equivalence can be given 
by the exact operator identity 



hvpk n cic n 

n=l 
flVp'K 



L 



u=i 

where the operators bi with I > 1 are defined as 



(1.3) 



h = Yj c Lc m +i (1.4) 



m=l 



and M = J2n= i c\c n is the fermionic particle number operator. The proof of 
the "Kronig identity" (1.3) is simple (see Ref. [4]) . The operators bi obey 



groundstate _•_ excited state 
20 

E E +20A F 



Figure 1.1. Example for the classification scheme for the excited states in terms of the numbers 
rij of upward shifts by j units of Aj?. In the example shown the nonzero rij are m = 1, n± = 
2, ns = 1 and m = 2. 



commutation relations [bi, by] = and for / > V 



(1.5) 



m=l 



For all iV-particle states |^ } ) = Jl^Li 4„l Vac ) in which the M (< N ) low " 
est one-particle levels are all occupied one obtains 



&,41iA m) > = <mA m) > 



(1.6) 



for < M, i.e. these operators obey boson commutation relations [bi, bj,] = 
5uA in this subspace of all possible N -particle states. 

Later it turns out to be useful to work with T = T — (T)o — n^N, where 
(T)o = Apnp(nF + l)/2 is the ground-state energy, fiQ = Af(uf + 1/2) 
is the chemical potential of the noninteracting fermions and M = M — ufI- 
Then f is of the form as the rhs of Eq. (1.3) with AT (AT + 1) replaced by M 2 . 
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2.2 Tomonaga (1950): Bloch's method of sound waves 
applied to interacting fermions 

When a two-body interaction between the fermions is switched on, the 
ground state is no longer the rilled Fermi sea but it has admixtures of (multi- 
ple) particle-hole pair excitations. In order to simplify the problem Tomonaga 
studied the high density limit where the range of the interaction is much larger 
than the interparticle distance, using periodic boundary conditions [5]. Then 
the Fourier transform v{k) of the two-body interaction is nonzero only for val- 
ues | A; | < k c where the cut-off k c is much smaller than the Fermi momentum 
k c <C kp. This implies that for not too strong interaction the ground state 
and low energy excited states have negligible admixtures of holes deep in the 
Fermi sea and particles with momenta \k\ — kp » k c . In the two intermediate 
regions around the two Fermi points ±Uf, with particle-hole pairs present, the 
dispersion is linearized in order to apply Bloch's "sound wave method" 

/c « ±/cp : Ek = Ep ±vp(k^-kp). (1.7) 

Tomonaga realized that the Fourier components of the operator of the density 

rL/2 . rL/2 

p n = / p(x)e~ %knX dx = \ ■ip\x)ip(x)e~ tknX dx 

J-L/2 J-L/2 
= E4 C «'+«> ( L8 ) 



where ct(c n ) creates (annihilates) a fermion in the state with momentum k n = 
^n, plays a central role in the interaction term, as well as the kinetic energy. 
Apart from an additional term linear in the particle number operator [4], which 
is usually neglected, the two-body interaction is given by 

V = 2l E v(kn)Pn P -n + ^ 2 5(0) (1.9) 

Tomonaga's important step was to split p n for \k n \ <C Hf into two parts, one 
containing operators of "right movers" i.e. involving fermions near the right 
Fermi point kp with velocity vf and "left movers" involving fermions near 
— k,F with velocity —vp 

Pn= 4i' C n'+n + ^ 4i' C n'+n = Pn,+ + Pn- (1-10) 

ri>0 n'<0 

where the details of the splitting for small \n'\ are irrelevant. Apart from the 
square root factor the p n a are similar to the hi defined in Eq. (1.4) . Their 
commutation relations in the low energy subspace are 

[Pm,a,Pn,/3] = OHnd a/3 S m - n i. (1.11) 
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If one defines the operators 



bn^^\^ + l° rn>0 n (1.12) 
1 1 1 p n for n < v 7 



and the corresponding adjoint operators b^ n this leads using pt a = p- n ,a to 
the bosonic commutation relations 

[b n ,b m }=0, [b n ,bl n ]=S mn l (1.13) 

Now the kinetic energy of the right movers as well as that of the left movers 
can be "bosonized" as in Eq. (1.3). The interaction V is bilinear in the p n as 
well as the p njQ . Therefore apart from an additional term containing particle 
number operators the Hamiltonian for the interacting fermions is a quadratic 
form in the boson operators 

n>0 ^ ^ ' 

, v( k n) fuiui i 7 , \1 , Tit; 



(bibi n + b_ n b n ) J + ^ [ vn n 2 + vjj 2 



2tt% 

= Hu + Htfj, (1.14) 

where N = N + + is the total particle number operator relative to the 
Fermi sea, J = J\f + — N- the "current operator", and the velocities are 
given by vn = vf + v(0)/ivh and v,j = vf- Here vn determines the en- 
ergy change for adding particles without generating bosons while vj enters the 
energy change when the difference in the number of right and left movers is 
changed. Similar to the discussion at the end of section 1.1 we have defined 
H = H — Eq — (po + v(0)no)AT, where Eq is the Hartree energy and uq the 
particle density. As the particle number operators N a commute with the boson 
operators 6 m (&m) the two terms Hb and Hj^ j in the Hamiltonian commute 
and can be treated separately. Because of the translational invariance the two- 
body interaction only couples the modes described by b\ L and b- n . With the 
Bogoliubov transformation a^ n = c n fejj — s n b- n the Hamiltonian Hb can be 
brought into the form 



Hb = ^2 Tv^nf^n a n + const., (1.15) 

where the uj n = vf \ k n \ y/l + vQc^JJttHvf follow from 2x2 eigenvalue prob- 
lems corresponding to the condition [Hb,^] = TkOnO^. For small k n one 
obtains for smooth potentials v(k) again a linear dispersion uj n ^ v c \k n \, with 
the "charge velocity" v c = ^/vnvj, which is larger than vf for v(0) > . 
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The expression for the coefficients c n and s n with (? n — s 2 n = 1 will be pre- 
sented later for the generalized model Eq. (1.17) . For fixed particle numbers 
N + and iV_ , the excitation energies of the interacting Fermi system are given 
by J2 m h^mUm with integer occupation numbers < n m < oo. For small 
enough excitation energies the only difference of the excitation spectrum for 
fixed particle numbers with respect to the noninteracting case is the replace- 
ment vf *-> v c . 

In his seminal paper Tomonaga did not realize the anomalous decay of cor- 
relation functions of the model because in his discussion of the density corre- 
lation function he missed the 2/cF-contribution discussed in section 3. 




-2.0 1 ■ 1 ■ 1 ■ 1 ■ 1 

-2.0 -1.0 0.0 1.0 2.0 

k/k F 

Figure 1.2. Energy dispersion as a function of momentum. The dashed curve shows the usual 
"nonrelativistic" dispersion and the full curve the linearized version used (apart from a constant 
shift) in Eq. (1.3) for k > for fixed boundary conditions. The dot-dashed parts are the addi- 
tional states for fco = — 1.5/cf- The model discussed by Luttinger corresponds to ko — > — oo. 



2.3 Luttinger (1963): no discontinuity at the Fermi surface 

Luttinger, apparently unaware of Tomonaga's work, treated spinless, mass- 
less fermions (in the relativistic sense, but c <-> vf) in one dimension, i.e. 
two infinite branches of right and left moving fermions with dispersion ±vpk 
[6]. As Luttinger himself made an error with the fact that his Hamiltonian 
is not bounded from below, it is useful to switch from Tomonaga's to Lut- 
tinger's model keeping a band cut-off ko such that k > ko = 2itmo/L with 
mo < for the right movers and correspondingly for the left movers (see Fig. 
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1.2). Fortunately Luttinger's error had no influence on his inquiry if a sharp 
Fermi surface exists in the exact ground state of the interacting model. After a 
rather complicated calculation using properties of so-called "Toeplitz determi- 
nants" Luttinger found that the average occupation (rik,+) in the ground state 
for k ^ kp behaves as 



where oil depends on the interaction strength (see below) [7]. "Thus, in this 
model, the smallest amount of interaction always destroys the discontinuity of 
(rifc) at the Fermi surface" [6]. This can be related to the fact that the equal time 
correlation functions (^ a (x)ip a (0)) decay as 1/| x\ 1+aL in the interacting sys- 
tem in contrast to (ip} i (x)i() a (0)} ~ (with d = 1) in the noninteracting 
case. Therefore oil is called the "anomalous dimension' '[8]. 

Apart from the different dispersion Luttinger also used a different interac- 
tion. In contrast to Tomonaga he only kept an interaction between the right and 
left movers but not the term ~ u(A; n )(6j l 6 ra + b^_ n b- n ) in Eq. (1.14) . In the 
limit of a delta interaction of the right and left movers his model is identical to 
the massless Thirring model (1958) [9] at that time not well known in the solid 
state physics community. 

2.4 Towards the "Luttinger liquid" concept 

Luttinger's treatment of the Dirac sea was corrected in a paper by Mattis 
and Lieb (1965) [10] which also offered a simpler way to calculate (nfc jQ ). 
The time dependent one-particle Green's function for the spinless Luttinger 
model was calculated by Theumann (1967) [11] by generalizing this method. 
She found power law behaviour in the corresponding spectral function p(k, lo), 
especially p(kp,uj) ~ aL\u\ aL ~ l , i.e. no sharp quasiparticle for k = kp 
consistent with Luttinger's result for the occupation numbers (Fig. 1.3). For a 
delta interaction her results agreed with an earlier calculation for the massless 
Thirring model by Johnson (1961) [12]. Later the time dependent one-particle 
Green's function was calculated by various other methods, e.g. using Ward 
identities (Dzylaloshinski and Larkin (1974) [13]) as well as the important 
method of the "bosonization of the field operator" (Luther and Peschel (1974) 
[14]) which will be addressed in detail in section 3. It was first proposed in a 
different context by Schotte and Schotte (1969) [15]. 

What is now usually called the "Tomonaga-Luttinger (TL) model" is the fol- 
lowing generalization of Eq. (1.14) 



k-k F aL 
k c 



sign(/c F - k), 



(1.16) 




n>0 



9i{kn) 

2nh 



) 
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+ 



92(k n . 

2irh 



(blbl n + b. n b n ) } + ^ {v N M 2 + vjJ 2 } , (1.17) 



where vn = vf + (54(0) + g2(0))/2irh and = vf + (54(0) - g2(0))/2irh. 
The interaction parameters g 2 (k n ) and gi(k n ) are allowed to be different. As 
Tomonaga's original model the TL model is exactly solvable, i.e. it can also 
be brought into the form of Eq. (1.15). The eigenvector components in cv n = 

Crfili ~ s nb- n are given by 



C n 





(1.18) 



withK n = y/vj(k n )/v N (k n ), where vj^ N) {k n ) = v F +[g i {k n )^g 2 {k n )]/2nh. 

The frequencies are given by oj n = \k n \^/vj{k n )vN{k n ) = \k n \v c {k n ). The 
TL-Hamiltonian corresponds to a fermionic Hamiltonian that conserves the 
number of right and left movers. 




-1 
(k-k F )/k c 



Figure 1.3. The full line shows the average occupation (n fcj+ ) for a TL model with ai = 0.6. 
The dashed line shows the expectation from Fermi liquid theory, where the discontinuity at 
determines the quasi-particle weight Zk F in p+(kF,uj). As discussed following Eq. (1.48) this 
can also be realized in a TL model with g-z (0) = 0. There also the details of the interaction are 
specified. 



A more general model including spin and terms changing right movers into 
left movers and vice versa is usually called the "g-ology model". An important 
step towards the general Luttinger liquid concept came from the renormaliza- 
tion group (RG) study of this model. It was shown that for repulsive interac- 
tions (see section 3) the renormalized interactions flow towards a fixed point 
Hamiltonian of the TL-type unless in lattice models for commensurate electron 
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fillings strong enough interactions (for the half filled Hubbard model discussed 
in section 4 this happens for arbitrarily small on-site Coulomb interaction U) 
destroy the metallic state by opening a Mott-Hubbard gap. The RG approach is 
described in detail in reviews by Selyom (1979) [16] and Shankar (1994) [17]. 
These results as well as insight from models which allow an exact solution 
by the Bethe ansatz led Haldane [18, 19] to propose the concept of Luttinger 
liquids (LL) as a replacement of Fermi liquid theory in one dimension, which 
"fails because of the infrared divergence of certain vertices it assumes to re- 
main finite" [19] . At least for spinless fermions Haldane was able to show that 
"the Bogoliubov transformation technique that solves the Luttinger model pro- 
vides a general method for resumming all infrared divergences present"[19]. 
Similar to Fermi liquid theory in higher dimensions this new LL phenomenol- 
ogy allows to describe the low energy physics in terms of a few constants, two 
for the spinless case: the "stiffness constant" K = Kq = ^vj/vn (also called 
g in various publications) and the "charge velocity " v c = ^Jvjvn- In his semi- 
nal paper Haldane showed explicitly that the LL relations survive in not exactly 
soluble generalizations of the TL model with a non-linear fermion dispersion. 
He also gave a clear presentation how to calculate general correlation functions 
and e.g. the occupancies shown in Fig. 1.3 for the TL model. The technical 
details are addressed in section 3. 

Before we do this two additional important aspects of LL-behaviour should 
be mentioned. The first concerns the strong influence of impurities on the low 
energy physics [20-25], especially the peculiar modification of the electronic 
properties of a LL when a single impurity with an arbitrarily weak backscat- 
tering potential is present. For a spinless LL with a repulsive two-body in- 
teraction, i.e. K < 1 a perturbative bosonic RG calculation [25] shows that 
the backscattering strength Vb is a relevant perturbation which grows as A^ 1 
when the flow parameter A is sent to zero. This leads to a breakdown of the per- 
turbative analysis in Vb- On the other hand a weak hopping between the open 
ends of two semi-infinite chains is irrelevant and scales to zero as 
Assuming that the open chain presents the only stable fixed point it was argued 
that at low energy scales even for a weak impurity physical observables behave 
as if the system is split into two semi-infinite chains. This leads to a conduc- 
tance which vanishes with a power law in T at low temperatures [25]. A more 
technical discussion is presented in section 3. 

Electrons are spin one-half particles and for their description it is necessary 
to include the spin degree of freedom in the model. For a fixed quantiza- 
tion axis the two spin states are denoted by a =|, [. The fermionic creation 
(annihilation) operators c\ l ± c {c nt ±^ a ) carry an additional spin label as well 
as the p n ,±,a and the boson operators b n ^ which in a straightforward way 
generalize Eq. (1.12). The interactions g v (k) with v = 2,4 in Eq. (1.17) 
become matrices g a v a in the spin labels. If they have the form g° a (k) = 
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5a,a'9i/\\(k) + 5a-a'9v±(k) it is useful to switch to new boson operators b n>a 
with a = c,s 

K,c = -j= (&»,T + Kl) 

K,s = ^(6„t-6 n ,|) , (1-19) 

which obey [6 a , n , & a , jn ,] = and b a>n , b\, n , = <5 aa /<5 n „/i. The kinetic energy 
can be expressed in terms of "charge" (c) and "spin" (s) boson operators using 
&n T^ ra 'T "nfinl = ^n,c^n,c + b\i,s^n,s- If one defines the interaction matrix 
elements g^ a {q) via 

5i/,c(?) = 9u\\(q) +9v±(q) 

9v, s {q) = 9v\\{q) - 9u±(q) , (1-20) 
and defines Af ajC (s) = (M*,! ^ ^a,i)/V^ one can write the TL-Hamiltonian 
H-Tif 1 f° r s P m one-half fermions as 

i^ /2) = H TL , C + Htl,s , (1.21) 

where the HrL,a are of the form Eq. (1.17) but the interaction matrix elements 
have the additional label a. The two terms on the rhs of Eq. (1.21) commute, 
i.e. the "charge" and "spin" excitation are completely independent. This is usu- 
ally called "spin-charge separation". The "diagonalization" of the two separate 
parts proceeds exactly as before and the low energy excitations are "massless 
bosons" uj nta m f a |fe n | with the charge velocity v c = {vj c vn c ) 1 I 2 and the 
spin velocity v s = {vj s vn s ) 1 ^ 2 ■ The corresponding two stiffness constants are 
given by K c = (v Jc /v Nc ) 1 / 2 and K s = {vjJv Ns ) l l 2 . Because of Eq. (1.21) 
the dependence of the velocities on the interaction strength (1.20) is obtained 
using the results for the spinless model following Eq. (1.18). 

The low temperature thermodynamic properties of the TL model including 
spin, Eqs. (1.17) and (1.21), can be expressed in terms of the four velocities 
vn c ,vj c , vn s ,vj s or the four quantities v c ,K c ,v s ,K s . Due to spin-charge 
separation the specific heat has two additive contributions of the same form as 
in Eqs. (1.1) and (1.2). If we denote, as usual, the proportionality factor in the 
linear T-term by 7 one obtains 

l = l(^ + V -Z), (1-22) 



7o 2 V v c v, 

where 70 is the value in the noninteracting limit. To calculate the spin suscep- 
tibility Xs one adds a term — hj\f s to B^jf 1 . Then by minimizing the ground 
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state energy with respect to N s one obtains (TV's) ~ h/vN s , i.e. Xs is inversely 
proportional to vn s - If one denotes the spin susceptibility of the noninteracting 
system by \sfi, this yields for the zero temperature susceptibility 

U ''' !<< — . (1.23) 



Xs,0 VN a Vs 

For spin rotational invariant systems one has K s = 1 [26]. The zero tem- 
perature compressibilty k is proportional to (d 2 Eo/dN 2 )~[ 1 which using Eqs. 
(1.17) and (1.21) leads to 

^ = ^ = K^. (1.24) 

v Nc v c 

A simple manifestation of spin-charge separation occurs in the time evo- 
lution of a localized perturbation of e.g. the the spin-up density. The time 
evolution a U:a (t) = a n ^ a e~ luJ "' at for a = c,s implies 

U (f \ — U f„2 -iu>„, a t _ 2 iu n , a t\ _}A ( '—iu„, a t _ iw n ,ai\ 
u n,a\^) — u n,a \^n,a e s ti,« e J "-n,a L n,ai>n,a \^ ) 

(1.25) 

If the initial state of the system involves a perturbation of right movers only, i.e. 
{b n ,a) = for n < and the perturbation is sufficiently smooth ((6 n ,a) / 
for < n <C n c only) the initial perturbation is split into four parts which 
move with velocities ±v c and ±v s without changing the initial shape. If only 
the initial expectation values of the b n j are different from zero one obtains for 
<5(p T (x,0)) = F(x) using Eq. (1.12) ' 



1 -±^F{X - Va t) + 1 —^F{X + v a t) 



(1.26) 



For the spin rotational invariant case K s = 1 there is no contribution which 
moves to the left with the spin velocity. Already for the pure 54-model with 
K c = 1 but v c / v s "spin-charge separation" of the distribution occurs. For 
the spinless model with §2 / the initial distribution splits into one right- and 
one left-moving part, which is often called "charge fractionalization" [27, 28]. 
Note that the splitting described in Eq. (1.26) is independent of the details 
of F(x) like the corresponding total charge. An additional comment should 
be made: spin-charge separation is often described as the fact that when an 
electron is injected into the system its spin and charge move independently 
with different velocities. This is very misleading as it is a collective effect of 
the total system which produces expectation values like in Eq. (1.26). 

The easiest way to understand the important manifestation of spin-charge 
separation in the momentum resolved one-particle spectral functions [29, 30] 
is to make use of the bosonization of the electronic field operators discussed in 
the next section. 
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3. Luttinger liquids - computational tools 

In section 2 many of the important features of LL's like the absence of a 
discontinuity at the Fermi surface were presented without giving any details 
how these properties are actually determined. As the most important tool the 
bosonization of the field operators is presented in detail in this section. This 
method is then used to calculate correlation functions like the one-particle 
Green's function and the 2 kp -density response function. In the second part 
of this section the TL model with additional interactions and (or) a one particle 
potential with a "backscattering" contribution is discussed. The model is no 
longer exactly solvable and RG arguments play an important role [16, 17, 25]. 

3.1 Bosonization of the field operator 

In the following a selfcontained presentation of the bosonization of a single 
fermion operator including a straightforward construction of the particle num- 
ber changing part ("Klein factor") is given. We present the bosonization of the 
field operator for the right movers described by the q + and just mention the 
corresponding result for the left movers. 

The starting point are the commutation relations the q + obey for m > 

[ftm,q,+] = — j=Q +mt+ , [£4,Q,+] = — ^Lq_ m>+ . (1.27) 



If (after taking the limit m — > — oo) one introduces the 27r-periodic auxiliary 
field operator ip+(v), where v later will be taken as 2ttx/L 

oo 

VM» EE Yl . d" 28 ) 

l=— oo 

it obeys the simple commutation relations 

[b m ,i> + (v)} = -_L e -<™V + (*;) ; [bl^ + (v)\ = -J^e imv ^ + (v) . 

v ?Yl \f 771 



(1.29) 

Products of exponentials of operators linear in the boson operators 

A + = J2Kbi ; B_ = J2fi n b n (1.30) 

with arbitrary constants \ n and /x„ obey similar commutation relations 

[b m ,e A +e B -] = \ m e A +e B - ; [bl, e A + e B -} = -fi m e A + e B ~ , (1.31) 
which follow from [b m , e Af> ™] = Ae Ab "\ We therefore make the ansatz 

i) + {v) = O+ivy^e^^, (1.32) 
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where the operator i<f)+(v) is given by [19] 

oo inv 

Then the operator 0+(v) commutes with all the b m and bl n . We next construct 
+ (v ) such that both sides of Eq. (1.32) yield identical matrix elements. 
As ip+(v) reduces the number of right movers by one, the operator + (v) 
also must have this property. In order to determine + (v) we work with the 
eigenstates of the noninteracting system ( the limit mo —> — oo is implied and 
np is an arbitrary positive integer later related to hp) 



\m,[ 



<n F +N- 

J 

-n- 

n=mo 

(1-34) 

It is easy to see that O + (v)\{0} b , N + , N-) has no overlap to excited states 



\{ mi } b ,N + ,N_)^H^=[ [I c[ 




{{m l } b ,N+ - \,N_\0+{v)\{Q} h ,N + ,N-) = 

(bi) m ' - ~ 
({0} b ,N + - l,N-\Yl^L=O + (v)\{0} b ,N + ,N.) . (1.35) 

As 0+(v) commutes with the b\ the rhs of Eq. (1.35) vanishes unless all m\ 
are zero. This implies 

6 + (v)\{0} b ,N + ,N^) =c + (N + ,N^,v)\{0} b ,N + -l,N_) , (1.36) 

where c + (N + , iV_ , v) is a c-number. In order to determine c + (JV+ , iV_ , v) we 
calculate ({0}{„ iV+ - 1, iV_|^ + (t>)|{0} 6 , 7V+, iV_) using Eq. (1.28) as well as 
Eq. (1.32). In the calculation of the matrix element with the fermionic form 
Eq. (1.28) we use Eq. (1.34) which yields 

({0} b ,N + - l,N-\ %+ \{0} b ,N+,N-) = (-lf-S ltnF+f f + . (1-37) 

The factor (— 1)^- occurs because we have to commute q i+ through the prod- 
uct of iV_ = — mo + \ + rip + iV_ fermionic operators of the left movers 
if we assume — mo + np to be odd. We note that no such factor occurs for 
the corresponding matrix element of the left movers. The calculation of the 
ground state to ground state matrix element of ip+(v) using Eq. (1.32) is sim- 
ple as both exponentials involving the boson operators can be replaced by the 
unit operator and the matrix element is just c + (N+, 7V_, v). The comparison 
therefore yields 

c + (N + ,N_,v) = (-l)^- e ™("F+iV + ) (L38) 
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and c_(iV+, N_, v) = e -iv(n F +N-\ Toge ther with Eq. (1.34) and the defini- 
tion M a = M a — (—mo + 1 + np)l this implies 

6 + (v)e-^ + ^(-lf-\{0} b ,N + ,N^) = \{0} b ,N + - 1,N_) . (1.39) 

If we apply the operator 6 + {v)e-^ nF+ ^ v [-1)^- to the states in Eq. (1.34) 
and use again that 0+(v) commutes with the boson operators we obtain 

d + (v)e- i ^ + ^(-lf-\{m l } b ,N + ,N^ = |{m,} 6 ,iV + - l,i\L> . 

(1.40) 

This shows that the operator U+ = + (v)e~ i{ - nF+ ^ v (-1)^- is independent 
of v and given by 

U + = £ Y>\{ m i}b,N+-l,N-){{m l } b ,N+,N-\. (1.41) 

N + ,N_ {mi} 

It follows immediately that U+ is unitary, i.e. U + Jj\_ = U\_U+ = 1. From Eq. 
(1.41) one can infer that for arbitrary functions / of the number operator M+ 
one has U+f(M+) = f(M+ + 1)E7+. 
To summarize we have shown that 

0+(v) = U + e^ nF+ ^+>{-\)^- . (1.42) 
In 6-{u) = [/_ e - i (™F+Ar_)« nQ factQr appears and therefore 6+(v) 

and 0_ (u) anticommute, which is necessary to enforce anticommutation rela- 
tions between ip + (v) and ijj- (u). It is an easy exercise to show that e.g. the an- 
ticommutation relations [ip + (v),tp + (u)] + = are fulfilled. In the calculation 
the properties of + (v) as well as the factor in Eq. (1.32) involving the boson 
operators enter. If one replaces the operators O a (v)e~ tav ^ a+nF ^ by "Ma- 
jorana fermions" r] a which commute with the boson operators and obey the 
anticommutation relations [r] a ,rjp\ + = 26 a pl, as often done in the literature, 
this yields [ip a (v),^ a (u)] + = [1 - cos (it - v )] e ^{u+v)(M a +n F ) ^ i e a viola . 
tion of the correct anticommutation relations. This implies that the U a have to 
be properly treated. In many publications they are written as U a = e 10 ™ , where 
the phase operators 9 a are assumed to obey the canonical commutation rela- 
tions (CCR) [M a , Oa] = *1 [19]. We do not use this concept here because no 
phase operator can be constructed which obeys the CCR as an operator identity 
[4,31-34]. 

In the following we will always use the "normal ordered" form (all boson 
annihilation operators to the right of the creation operators) of the bosonization 
formula Eqs. (1.32, 1.33). Alternatively one introduces a convergence factor 
e~ nX / 2 , whith A —> and works with the Hermitian Bose fields $ a (v) = 
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<t>a(v) + ^a(v) as well as the fields $ + ± The derivatives of the latter 
fields are related to the total current and the deviation of the total charge density 
from its average value [35]. As we work with an interaction cut-off k c the 
introduction of A is not necessary and because of the space limitation this field- 
theoretic formulation is not used here. 

3.2 Calculation of correlation functions for the TL model 

In order to calculate correlation functions of the TL model with nonzero in- 
teractions it is necessary to express the field operator V>+(t>) Eq. (1.32) in 
terms of the a„,o4 instead of the 6 n ,&t because the former have a simple 
time dependence and for the temperature dependent expectation values one 
has (a) m a n ) = 8 mn riB(oj n ), where nsioj) = l/(e l3uJ — 1) is the Bose func- 
tion. For the ground state calculation all one needs is a n |$o) = without 
using the explicit form of the interacting ground state |<I>o)- 

For periodic boundary conditions one has b m = c m a m + s m a^ m where 

the operators a m and al m commute. Therefore e^+^ (and e l< ^+^) in Eq. 
(1.32) can be written as a product of two exponentials with the annihilation 
operators to the right. After once using the Baker-Hausdorff formula, e A+B = 
e A e B e~^ A,B ^ if the operators A and B commute with [A, B], in order to com- 
plete the process of normal ordering one obtains for the physical field operator 
V'a (x) = 4> a (27rx /V)/ \fZ for a system of finite length L with periodic bound- 
ary conditions [36] 

^ Q ( x ) = ^-6 a {^pj e *!.(*) e ixa(*) (1.43) 

with 

iXa{x) = ^7j=r (c m e ikmX a m - s m e- ikmX a. m ) , (1.44) 

. . , , _ n=i Sn/n and 9{x) is the unit step function. 

This is a very useful formula for the calculation of properties of one-dimensional 
interacting fermions. For the special choice s n = s(0)e~ n ^ nc where n c = 
k c L/2ir is determined by the interaction cut-off, A(L) can be calculated ana- 
lytically using Y^?=i z n /n = — log (1 — z). For n c » 1 this yields A(L) = 
(4ir/k c L) s (°) which shows that the prefactor in Eq.(1.44) has an anomalous 
power law proportional to (1/L)5+ S (°). This implies that the c n a scale like 
(1/L) s2 (°). 

The time dependent operator i/j + (x,t) follows from Eq. (1.44) by replacing 

a m and a_ m by a m e~ luJmt and a.- m e~ tu)mt and U + in 0+ by U+(t). Various 
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kinds of time dependent correlation functions can quite simply be calculated 
using this result. Here we begin with iG+(x, t) = (^+(0, 0)tp + (x, t)}. 
As U + commutes with the bosonic operator the particle number changing op- 
erators lead to a simple time dependent factor 

ulu+(t)\§ (N + ,N-)) = e - i i Eo ^+^- ) - Eo ^+' 1 ^- ) V\^ (N + ,N^)). 

(1.45) 

As tp+( x ) m Eq.(1.43) is normal ordered in the a's one has to use the Baker- 
Hausdorff formula only once to normal order ip\_(Q,Q)ip + (x,t). This yields 
with hp = 2irnp/L 

ie l ^G<{x,t) = ^!i^ e ^ e [x+ (o,o)4 (*,*)] (1.46) 
L 

L 

where \i = E (N+,N-) - E (N+ - 1,7Y_) is the chemical potential. The 
analytical evaluation of the sum (integral in the limit L — ► oo) in the exponent 
in Eq. (1.46) is not possible. An approximation which gives the correct large 
x and t behaviour [37] is to replace oo n by v c k n . This yields for L — > co with 
the exponential cut-off for the s n used earlier [14] 

• ifit < -">■ eikpX ' r<2 ' S2(0) 

+ ^ ' ^ 2ir x — v c t — iO (x — v c t — ir)(x + v c t + ir) 

(1.47) 

where r = 2/k c . As (^.(0,0)ip + (x,0)) for large x decays proportional to 
(l/x) 1+2s (°) the anomalous dimension for the spinless model is given by 

a L = 2s 2 (0) = (K — l) 2 /2K. (1.48) 

Luttinger's result for (nk,+) follows by performing the Fourier transform with 
respect to x. The full line in Fig. 1.3 was calculated with s 2 n = 0.3e~ 2fcn / fcc , 
while the dashed curve corresponds to s 2 = 0.6(k n /k c )e~ 2kn ^ kc . The latter 
example corresponds to an interaction with §2(k — > 0) —>■ which leads to a 
vanishing anomalous dimension aL- In this case the occupancies have 
a discontinuity at kp as in a Fermi liquid [38]. An efficient numerical algo- 
rithm to calculate (rik,+) for arbitrary forms of s 2 t is described in the appendix 
of reference [2]. 

The spectral function p < (k, oo) relevant for describing angular resolved pho- 
toemission is obtained from Eq. (1.47) by a double Fourier transform 

p<{k,uj) = (cl + S[oj + (H-E (N + -l,N-))]c k , + ) (1.49) 



1 roo POO 

= — dte iujt / dxe- ikx ie^G<{x,t). 

■^TT J — oo J — oo 
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As Eq. (1.47) is reliable in the large x and t limit its use in Eq. (1.49) correctly 
describes the spectral function for k m kp and u> « [39]. Using the variable 
substitutions u T = x^fv c t the double integral factorizes and with the additional 
approximation iO — > ir on the rhs of Eq. (1.47) one obtains [11, 14] 

p<(k F + k,uj) ~ e{-u-v c \^)(-u+kv c )^- l {-u-lv c )^r e r "/ v * . (1.50) 

Without the additional approximation there is an additional weak dependence 
on u> + kv c [29]. The complete spectral function p + {k,uj) = p+{k,uj) + 
p^{k,uj), where p^(k,uj) is defined via iG+(x,t) = (ip + (x,t)ip\_(0,0)) can 
be obtained using p+(/cf + k,uj) = /?+(&f — k, —uj) which follows from the 
particle-hole symmetry of the model. The absence of a sharp quasi-particle 
peak is manifest from p + (kp,tu) ~ ap\uj\ aL ~ 1 e~ r \ u> ^ Vc . 

In order to calculate correlation functions of the spin one-half TL model 
the operators 6 n (7 which appear in the generalization of the bosonization for- 
mula Eqs. (1.32) and (1.33) have to be replaced by the spin and charge bosons 
b n ,a = (b n ,c + &bn tS )/\/2. Because of the exponential occurence of the bo- 
son operators in Eq. (1.32) and spin-charge separation Eq.(1.21) the Green's 
function a (x, t) factorizes into a spin and a charge part, which both are 
of the form as the square root of the function on the rhs of Eq. (1.47) . This 
square root results from the factors l/\/2 in the expression for the 6„ i(7 . In 
the spin factor the charge velocity v c is replaced by the spin velocity v s . For 
the average occupation numbers one again obtains Luttinger's result Eq. (1.16) 
with ai = «c(0) + = «c + « s - The individual contributions can be 

expressed in terms of the K a = {v,j a /v Na ) 1 / 2 as a a = (K a - l) 2 /(AK a ). 

As in the spinless model the fermionic (creation) annihilation operators Cn} a ,a 
scale like {l/L) aL / 2 . For spin rotational invariant systems one has K s = 1, 
i.e. no contribution to the anomalous dimension cil from the spin part of the 
Hamiltonian [26]. For the momentum integrated spectral functions one ob- 
tains p ata (uj) ~ \u\ aL as in the spinless model [40]. The /c-resolved spectral 
functions p a ,a{k, on the other hand show a drastic difference to the model 
without spin. The delta peaks of the noninteracting model are broadened into 
one power law threshold Eq. (1.47) in the model without spin and two power 
law singularities (see Fig. 1.4) in the model including spin [29, 30, 37] (for 
ai < 1/2 in the case of a spin independent interaction). The "peaks" disperse 
linearly with k — kp- 

It is also straightforward to calculate various response functions for the TL 
model. We discuss the density response function R(q, z) = —({p q ; f>- q )) z / L 
of the spinless model for q m and q rs ±2kp, where 



((A;B)) z = --j o ([A(t),B]y zt dt 



(1.51) 
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Figure 1.4. Spectral function p +tt7 (kF + k,ui) as a function of normalized frequency for 
k — —k c /10 for the TL- model with a spin independent interaction. The parameters are chosen 
such that v c = 2vf and o?l=1/8. 



involves the retarded commutator [41] and z is a frequency above the real axis. 
From the decomposition [42] i/)(x) ~ 4>+{x) + ip-(x) of the field operator 
ip{x) in the original Tomonaga model it is obvious that the operator p(x) = 
ip^(x)ip(x) of the density (see Eq. (1.8) ) has two very different contributions 

p(x) « p + (x) + p-(x) + (i)\(x)il)-(x) + /i.e.) (1.52) 

= Po( X ) +P2k F (x). 

The spatial Fourier transform of p is linear in the boson operators Eq. (1.12) 
and the q w contribution of the density response function [R(q, z)]q defined 
with the operators (po)q follows using the (linear) equations of motion for the 

b n (t) as 

[% * = ^Ki^ (1 - 53) 

This exact result for the q m contribution agrees with the RPA result for 
the original Tomonaga model. This fact, not mentioned in Tomonaga's paper 
[5] as the RPA paper by Bohm and Pines [43] was not yet published, was 
"discovered" many times in the literature. For the spin 1/2-model [R(q, z)]q 
has an additional factor 2 and one has to replace v j by vj c . 

The real part of the (q = 0) frequency dependent conductivity a (to + iO) 
follows from [R(q, u> + zO)]o by multiplication with ie 2 uj/q 2 and taking the 
limit q — ► 0. This yields for the spinless model 



(h/e 2 )Rea(oj + iO) = v j5{uj) = Kv c 5{uj) (1.54) 
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For the Galilei invariant Tomonaga model Eq. (1.14) one has vj = vf, i.e. the 
weight D of the zero frequency "Drude peak" is independent of the interaction, 
as expected. As D apart from a constant is given by the second derivative of 
Eq($)/L with respect to the magnetic flux through the Id ring [44], K (or K c ) 
can be obtained from a ground state calculation for microscopic lattice mod- 
els using Ki c \ = {Dn/ DqKq) 1 / 2 , where k is the compressibility discussed in 
Eq. (1.24). The anomalous decay of the correlation functions for these mod- 
els, which are more difficult to calculate directly, can then be quantitatively 
predicted if Haldane's LL concept is taken for granted. For a weak two-body 
interaction the result for Kt c \ — 1 linear in the interaction follows from first 
order perturbation theory for the ground-state energy, which involves the (non- 
selfconsistent) Hartree and Fock terms. As they are independent of the mag- 
netic flux, D/Dq has no term linear in v, i.e. K c m (k/kq) 1 / 2 = (vf/vNc) 1 ^ 2 , 
which holds exactly for Galilei invariant continuum models [45]. Performing 
the second derivative of (N) with respect to N yields [46] 

zirhVF 

In the spinless case the factor 2 in front of v(0) is missing in the result for 
K. Instead of D as the second input besides k one can obtain v c directly by 
calculating the lowest charge excitation energy (see section 4). 

The easiest way to calculate the q ±2kp contribution to the density re- 
sponse is to use the bosonization of the field operators [14]. The first step is to 

normal order ip' + (x)ip-(x) using Eq. (1.43) . This gives a factor e^- x+ ^' x -^ 

which using [x+(a;), X-( x )] = ~~ 2 J2 m >o c-mSm/m together with the factor 
A 2 (L) leads to 

*Xix)*-{x) = ^ (^)^t(^)0-(^) e -^ t ^ e -^ ) (1.56) 



with 



Ax(» = X+O) - X-( x ) = ~ 




Here ao is a dimensionless constant of order unity and the exponent K — 1 
of the second factor on the rhs follows using 2s^ + 2c m s m = K m — 1. The 
importance of this factor for impurity scattering in Luttinger liquids was first 
pointed out by Mattis (1974) [21] and will be discussed later. The calculation 
of the two terms of the commutator ([ip' + (x,i)ip-(x,i), ip_(0, 0)^ + (0, 0)]) is 
then straightforward and one obtains for the spectral function of the q ps ±2fcp 
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response function the power law behaviour [14] 

~ ~2fc2 ) (1-57) 

The static ±.2Uf + <2 response diverges proportional to |Q| 2 ^ -1 ^ which has 
to be contrasted with the logarithmic singularity in the noninteracting case. In 
the model including spin the exponent 2K — 2 is replaced by K c + K s — 2. 

The pair propagator P(q, uj) resulting from the response function for A = 
ip\_(x)ip^_(x) and B = ip-(0)ip+(Q) was found by Luther and Peschel to be 
the same as the 2/cp-density response, provided the sign of the interaction is 
reversed [14]. An attractive interaction leads to a power law divergence in 
P(q = 0, lo = 0) as the temperature is lowered, indicative of large pairing 
fluctuations. 



3.3 The TL model with additional interactions and 
perturbations 

The exact solution of the TL model essentially depends on the fact that the 
numbers of right and left movers are conserved. This symmetry can be de- 
stroyed by a one-particle potential with ±2/cF-Fourier components or by inter- 
action terms which change the individual particle numbers, like 2fcp-"backscattering" 
terms or Umklapp-terms for a half-filled band. With such additional terms the 
model is in general no longer exactly solvable. Important insights about the in- 
fluence of such terms have come from a (perturbational) RG analysis [16, 25]. 

3.3.1 Impurity in a spinless TL model 

We begin with the spinless model with an additional impurity which is de- 
scribed by 

Vi = J [V F (x)po(x) + V B {x)p 2kF {x)] dx = V F + V B , (1-58) 

where Vp describes the forward and Vb the backward scattering due to the im- 
purity and the two different operators for the densities are defined in Eq. (1.52) 
. As the forward scattering term is linear in the boson operators it can be treated 
in an exact way. The backscattering term has the property [Vb,No\ / and 
the model can no longer be solved exactly (except for K = 1/2 and a spe- 
cial assumption about Vb, as discussed below). For a zero range impurity it 
follows directly from Eq. (1.56) that Vb scales as (1/L) K while Htl in Eq. 
(1.17) scales as 1/L. Therefore the influence of V B depends crucially on the 
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sign of the two-body interaction [21, 22]. For repulsive interactions one has 
K < 1 which shows that Vb is a relevant perturbation. For K > 1, i.e. an 
attractive interaction, is irrelevant. A detailed RG analysis of the problem 
was presented in a seminal paper by Kane and Fisher [25]. For a zero range 
backscattering potential and two-body interaction they mapped the problem to 
a local bosonic sine-Gordon model [25, 35, 47]. The subsequent RG analysis 
shows that the backscattering amplitude scales as A^" -1 when the flow param- 
eter A is sent to zero [48], as can be anticipated from Eq. (1.56) . This leads to 
the breakdown of the perturbational analysis in Vb for repulsive interactions. 
As already mentioned in section 2 this analysis was supplemented by a RG 
analysis of a weak hopping between two semi-infinite chains. The weak hop- 
ping scales to zero like A aB for repulsive interactions, where as = K^ 1 — 1 is 
the boundary exponent. It describes e.g. the different scaling p(x, lo) ~ \u\ aB 
of the local spectral function near a hard wall boundary of a LL [25, 49, 50]. 
These scaling results together with the asumption mentioned in section 2 leads 
to the "split chain scenario" in which even for a weak impurity the observ- 
ables at low energies behave as if the system is split into two chains with fixed 
boundaries at the ends. Within the bosonic field theory this assumption was 
verified by quantum Monte Carlo calculations [51] and the thermodynamic 
Bethe ansatz [52]. 

This implies e.g. for the local density of states p(x,uj) ~ \uj\ aB for small 
and x near the impurity like in a LL near a hard wall. The transmission 
through the impurity vanishes near Uf proportional to ~ \k — kp\ 2oLB which 
leads to a conductance G{T) which vanishes with temperature T in power law 
fashion G(T) ~ T 2aB [25]. 

Additional insight comes from the analysis for the special value K = 1/2 
[25, 35, 32] . For V B {x) = V B S(x) the expression for Ax(0) in Eq. (1.56) can 
be written in terms of new boson operators a m = (a m — a_ m ) / \/2. If one ne- 
glects the momentum dependence of K m in Eq. (1.56) and puts K m = 1/2 
one obtains iAx(O) = Em>i ttm/\/™ as in the bosonization of a single 
field operator Eqs. (1.32) and (1.33) . It is then possible to refermionize the 
K = 1/2-TL model with a zero range impurity. Even the Klein factors can 
properly be handled [32] and one obtains a model of "shifted noninteracting 
Fermi oscillators" which can be solved exactly introducing an auxiliary Ma- 
jorana fermion [35, 32]. Unfortunately the local densities of states cannot be 
calculated exactly because of the complicated nonlinear relationship between 
the original fermion operators and the fermion operators which diagonalize 
the shifted Fermi oscillator problem [32]. Additional results for the transport 
through a spinless LL containing one impurity were obtained by mapping the 
problem onto the boundary sine-Gordon model and using its integrability [53]. 

In order to bridge the two regimes treated by Kane and Fisher one can use a 
fermionic RG description bearing in mind that it is perturbational in the two- 



23 



body interaction [54, 55]. It shows that the long range oscillatory effective 
impurity potential is responsible for the "splitting", for site impurities as well 
as for hopping impurities of arbitrary strength. For realistic parameters very 
large systems are needed to reach the asymptotic open chain regime [55]. 
Hence only special mesoscopic systems, such as very long carbon nanotubes, 
are suitable for experimentally observing the impurity induced open boundary 
physics. 

For a discussion of the impurity problem in the TL model including spin see 
also reference [56]. 

3.3.2 The TL- model with additional two-body interactions 

Tomonaga was well aware of the limitations of his approach for more generic 
two-body interactions ("In the case of force of too short range this method 
fails"[5]). We therefore first discuss Tomonaga's continuum model in this short 
range limit k c » hp opposite to the one considered in section 2. Then low en- 
ergy scattering processes with momentum transfer « ±2kp are possible and 
have to be included in the theoretical description of the low energy physics. 

In the "g-ology" approach one linearizes the nonrelativistic dispersion around 
the two Fermi points and goes over to right- and left-movers as in section 2. 
Then the "2fcir"-processes are described by the additional interaction term 



x. 



<J,cr' 

(1.59) 

For a spin-independent two particle interaction one has gn = g±± = g\. For 
the zero range interaction assumed in Eq. (1.59) one has to introduce band cut- 
offs to regularize the interaction term. The RG flow equations for the cut-off 
dependent interactions on the one-loop level are quite simple [16]. If s runs 
from zero to infinity in the process of integrating out degrees of freedom one 
obtains for spin-independent interactions 



dgi (f) 1 

ds itHvf 
dg 2 {s) 1 



gi(s) (1.60) 
gi(s) 



ds 2tt%vf 

and g 4 is not renormalized. Obviously g\ (s) can be obtained from the first 
equation only 

1 b whv F 

where gi is the starting value. It is easy to see that gi(s) — 2g 2 (s) = gi — 2g 2 
holds by subtracting twice the second equation from the first in Eq. (1.60) . 
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In the following we use the notation g* = g u {s — > oo). Now one has to 
distinguish two cases: 

for gi > one renormalizes to the fixed line g\ = 0, g\ = g2 — gi/2 and 
the fixed point Hamiltonian is a TL model which shows the generic importance 
of the TL model for repulsive interactions. In this case the g\ -interaction is 
called marginally irrelevant. For the nonrelativistic continuum model with a 
spin independent interaction one has g\ c = 2v(0) — v(2kp) and g^ s = and 
for the stiffness constant K c = [(2-kvf + g\ c - #| c ) / (2ttv f + g\ c + 92c)} X ^ ~ 
1 — [2u(0) — v(2kp)]/(27rhvF) and K s = 1. Due to the approximations made, 
also here only the result for K c — 1 linear in v is reliable. The agreement 
with the direct calculation Eq. (1.55) shows explicitly to leading order in the 
interaction that Haldane's Luttinger liquid concept is consistent. 

For g\ < the solution (1.61) shows that gi(s) diverges at a finite value 
of s. Long before reaching this point the perturbational analysis breaks down 
and all one can say is that the flow is towards strong coupling. In this case the 
<7i- interaction is called marginally relevant. In order to to obtain an under- 
standing of the strong coupling regime it is useful to bosonize the additional 
interaction in Eq. (1.59) [57]. The term proportional to gn is of the form 
of a g 2 1 1 -interaction and therefore bilinear in the boson operators Eq. (1.19) . 
For the gi^-term one uses the bosonization of the field operators Eqs. (1.32) 
and (1.33) with additional spin labels. As the g^-term contains field opera- 
tors ip^(x)ip a i{x) of opposite spin it only involves "spin bosons" Eq. (1.19) 
, which implies "spin-charge separation" also for this model [58]. The charge 
part stays trivial with massless charge bosons as the elementary interactions. 
Luther and Emery showed that for a particular value of g-y\\ the gi^-term can 
be written as a product of spinless fermion field operators and the exact so- 
lution for the spin part of the Hamiltonian is possible using refermionization 
[57], discussed earlier in connection with the backscattering impurity. The 
diagonalization of the resulting problem of noninteracting fermions is simple 
and shows that the spectrum for the spin excitations is gapped. It is generally 
believed that these properties of Luther-Emery phases are not restricted to the 
solvable parameter values. 

Strong coupling phenomena which lead to deviations from LL-properties 
with gapped phases are discussed in detail in section 4 for lattice models. There 
in case of commensurate filling Umklapp processes can become important, e.g. 
for half filling where two left movers from the vicinity of the left Fermi point 
are scattered into two right movers near the right Fermi point or vice versa. As 
G = 4:kp is a reciprocal lattice vector such a scattering process is a low energy 
process conserving quasi-momentum. In the g-ology model such processes are 
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described by an additional term 





,2ik F (x+y) 



+H.c] dxdy 



(1.62) 



Umklapp processes for a = a' are only possible for nonzero interaction range. 

4. Results for integrable lattice models 

As mentioned in subsection 2.4, results for integrable models which can 
be solved exactly by the Bethe ansatz played a central role in the emergence 
of the general "Luttinger liquid" concept [18]. It is therefore appropriate to 
shortly present results for the two most important lattice models of this type, 
the model of spinless fermions with nearest neighbour interaction and the ld- 
Hubbard model. (We put fi = 1 in this section.) 

4.1 Spinless fermions with nearest neighbour interaction 

The one-dimensional single band lattice model of spinless fermions with 
nearest neighbour hopping matrix element t(> 0), and nearest neighbour in- 
teraction U (often called V in the literature) is given by 



where j denotes the sites and the hj = c -Cj are the local occupation number 
operators. In the noninteracting limit U = one obtains for lattice constant 
a = 1 the well known dispersion e& = — 2t cos k. For the following discussion 
of the interacting model (U / 0) we mainly focus on the half filled band case 
hp = vr/2 with vf = 21 In contrast to the (continuum) Tomonaga model 
Umklapp terms appear when the interaction term in Eq. (1.63) is written in the 
^-representation. As discussed below they are irrelevant at the noninteracting 
(U = 0) fixed point [17]. Therefore the system is a Luttinger liquid for small 
enough values of \U\. The large U limit of the model is easy to understand: 
For U » t charge density wave (CDW) order develops in which only every 
other site is occupied thereby avoiding the "Coulomb penalty". For large but 
negative U the fermions want to be as close as possible and phase separation 
occurs. For the quantitative analysis it is useful that the model in Eq. (1.63) can 
be exactly mapped to a S = 1/2-Heisenberg chain with uniaxially anisotropic 
nearest neighbour exchange ("XXZ" model) in a magnetic field by use of the 
Jordan- Wigner transformation [59]. For U > this model is also called the an- 
tiferromagnetic Heisenberg-Ising model. The point U = U c = 2t corresponds 
to the isotropic Heisenberg model. For U > 2t the Ising term dominates and 



'3+1 



+ H.c. 



■)+Uj2 njh j+ i = f + U, (1.63) 
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the ground state is a well defined doublet separated by a gap from the contin- 
uum and long range antiferromagnetic order exists. For — 2t < U < 2t there 
is no long range magnetic order and the spin-excitation spectrum is a gapless 
continuum. The mapping to the XXZ-model therefore suggests that the spin- 
lesss fermion model Eq. (1.63) in the half filled band case is a Luttinger liquid 
for |Z7| < 2t. 

Before we present the exact results for the Luttinger liquid parameters K 
and v c from the Bethe ansatz solution [18, 60], we shortly discuss the RG 
approach to the model. A perturbative RG calculation around the free fermion 
fixed point is discussed in detail in Shankar's review article [17]. The first step 
is to write the four fermion matrix elements of the interaction U in Eq. (1.63) in 
the k -representation. This yields for a chain of N sites with periodic boundary 
condition and values kj = 2-irj/N in the first Brillouin zone 

(k 1 ,k 2 \U\k 3 ,k 4 ) = 2Ucos (h ^) ^2 d kl+ k 2 ,k 3 +k4+277m (1-64) 

m=0,±l 

The m = term on the rhs of Eq. (1.64) represents the direct scattering terms 
and the m = ±1 terms the Umklapp processes. The matrix element anti- 
symmtrized in ks and k± is proportional to sin \{k\ — &2)/2] sin [(kz — k 4 )/2\. 
Therefore the low energy Umklapp Hamiltonian scales like (1/L) 3 which shows 
that it is strongly irrelevant at the free field fixed point [17]. This analysis 
confirms the Luttinger liquid behaviour for small values of U, but gives no 
hint about the critical value U c for the CDW transition. With the separation 
U = Uq + ^umklapp implied by Eq. (1.64) one can do better by first treating 
T + Uq by bosonization and then perform the RG analysis around the corre- 
sponding TL fixed point to get information for which value of U the Umk- 
lapp term starts to be a relevant perturbation. For this analysis it is easier 
to work directly with the unsymmetrized matrix elements in Eq. (1.64). As 
k\ — &3 ~ ±7r for the low energy Umklapp processes this leads after extending 
the (linearized) dispersion of the right and left movers from —00 to 00 to a 
53 -interaction with a range of order r = a. The scaling dimension of the cor- 
responding H-j2 follows using bosonic normal ordering as in Eq. (1.56). For 
x — y of order r or smaller one obtains 

^ + {x)^ + {y)^{y)^(x)L 2 ~ (^p)' {j)^ ( L65 ) 

x (U\_ ) 2 Ut e 2ikF ( x+y) e iBi (:E ' y) e iB{x ^ , 

where B(x,y) = X-( x ) + X-(v) ~ X+( x ) ~ X+(v) with Xa(x) defined in 
Eq. (1.44). The first factor on the rhs is due to the Pauli principle and describes 
the same physics as the two sin-factors mentioned above for small arguments. 
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Therefore the second factor has to provide more than two powers of L to make 
the Umklapp term a relevant perturbation, which happens for K < 1/2. As 
discussed below, the exact Bethe ansatz result for K yields U c = 2t. If one 
uses the simple linear approximation for K — 1 in Eq. (1.55) one obtains with 
Eq. (1.64) K hn = 1 - U/(irt) for the critical value U^ a /t = vr/2, not too bad 
an approximation. 

Exact analytical results for the Luttinger liquid parameters for the half filled 
model can be obtained from the Bethe Ansatz solution [18, 60, 61]. It is not 
necessarary to address the anomalous decay of the correlation functions di- 
rectly, but one can use a ground state property and the lowest charge excitation 
to extract the parameters, as was dicussed in connection with Eq. (1.55). This 
yields for the stiffness constant K = ir/{2 arccos {—U/2t)\ and for the charge 
velocity v c = — (U/2t) 2 /[it — arccos (—U/2t)]. For repulsive interac- 

tions U > the value of K decreases monotonously from the noninteracting 
value K = 1 to K = 1/2 for U = 2t, which corresponds to an anomalous di- 
mension o>l = {K+\/K)/2 — \ = 1/4. For attractive interactions K diverges 
when U approaches —2t, and the charge velocity v c goes to zero. Results for 
the Luttinger liquid parameter K for less than half filled bands are shown in 
Fig. 1.5 [62]. The limit a — > and n — > corresponds to the continuum limit. 




Figure 1.5. Luttinger liquid parameter K from the Bethe ansatz solution as a function of the 
band filling n for different values of U (t = 1). The short dashed curve shows the infinite U 
result (1/2 + |n-l/2|) 2 . 



As the interaction goes over to a contact interaction its effect vanishes because 
of the Pauli principle and K goes to 1. For small enough values of U the linear 
approximation Eq. (1.55) K hn = 1 — ?7sin (nrf/irt provides a good approx- 
imation for all values of n, in contrast to the Hubbard model discussed below. 
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In the infinite U limit the Bethe Ansatz equations simplify considerably and 
the ground-state energy as well as low lying excited states can be calculated 
analytically [61]. With these results it is easy to show that K = (1 — n) 2 holds 
for < n < 1/2, i.e. K = 1/4, is the lower bound for K in the LL regime of 
the model [18]. The corresponding upper bound of the anomalous dimension 
is = 9/8. In order to achieve larger values of oil the model in Eq. (1.63) 
has to be generalized to include longer range interactions [63]. 

4.2 The Hubbard model 

As there exists an excellent review on the LL behaviour in the ld-Hubbard 
model [64], the following discussion will be rather short. As the model in- 
cludes spin the on-site interaction between electrons of opposite spins is not 
forbidden by the Pauli principle. This is taken as the only interaction in the 
model. The Id Hubbard Hamiltonian reads 

H = ~ t J2 ( c Ls+i, CT + H - c ) ■ r Y, "j. ".!..- d-66) 

In the extended Hubbard model a next nearest interaction term VJ2j ^j^j+i 
with hj = hjj + hj : i is added [65]. In order to show the important difference 
to the spinless model Eq. (1.63) we again first discuss the half-filled band case, 
which is metallic for U = 0. For U S> t the "Coulomb penalty" is avoided 
when each site is singly occpied. Then only the spin degrees of freedom matter. 
In this limit the Hubbard model can be mapped to a spin- 1/2 Heisenberg anti- 
ferromagnet with an exchange coupling J = At 2 /U. In the charge sector there 
is a large gap A c ~ U while the spin excitations are gapless. The Id Hubbard 
model can also be solved exactly using the Bethe ansatz [66] and properties 
like the charge gap or the ground-state energy can be obtained by solving Lieb 
and Wu's integral equation. In contrast to the spinless model described in 
the previous subsection the charge gap in the Hubbard model is finite for all 
U > 0. While for U 3> t it is asymptotically given by U it is exponentially 
small, A c « (8t/ir)^/TT/texp (-2nt/U), for < U < t. This shows that 
the Umklapp term is no longer irrelevant at the free field fixed point. The Pauli 
principle factor of Eq. (1.65) is missing as the interaction is between electrons 
of opposite spin. The Umklapp term is therefore a marginal perturbation. The 
RG analysis [16] shows that the Umklapp term is marginally relevant while 
the 2fcir-backscattering ("51") interaction is marginally irrelevant for U > as 
discussed following Eq. (1.60). 

When the band is not half filled Umklapp is not a low energy process and 
the Hubbard model is a Luttinger liquid with K s = 1. The LL parameters K c 
and v a can be obtained by (numerically) solving Lieb an Wu's integral equation 
[67]. Even for < U <C t the perturbative result Eq. (1.55) works well only for 
intermediate filling n = N e i/N w 0.5, where N e i is the number of electrons 
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(half filling corresponds to n = 1) . In the limit n — > the Fermi velocity 
■up = 2t sin (irn/2) goes to zero but 2u(0) — v(2kp) = U stays finite and the 
correction term increases with decreasing n in contrast to the spinless model. 
The Bethe ansatz results show that K c — ► 1/2 for n — > as well as n — ► 1 
for a/Z [/ > 0. For U — > oo it leads to i<C c — ► 1/2 for all fillings n different 
from 1. In this limit the velocities are given by v c = 2tsin(7rn) and v s = 
(27rt 2 /U)[l - sin(27rn)/(27rn)], i.e. the spin velocity goes to zero [64, 67]. 
The U = oo results for v c and K c can be understood without the Bethe ansatz 
solution. Double occupancies of the lattice sites are forbidden and the system 
behaves like a system of noninteracting spinless fermions with fcp replaced by 
2kF [64]. The spin degrees of freedom play no role and any configuration of 
the spins gives an eigenfunction of the same energy. This immediately explains 
the result for v c mentioned above. For a TL model with spin one obtains (for 
fixed iV T - N0 fromEqs. (1.17) and (1.21) L(d 2 E /dN 2 ) L = nv Nc /2, while 
the factor 1/2 is missing in the spinless case. The formula for the spinless 
case can be used to calculate L(d 2 Eo/dN 2 )L for U = oo with vn replaced 
by VF{2kp), using the spinless fermion analogy. This yields vn c = 2v c i.e. 
K c = l/2. 

As the calculation of correlation functions not only requires excitation en- 
ergies but also many electron matrix elements which are difficult to evaluate 
using the Bethe ansatz, various numerical methods have been used to study 
e.g. the manifestation of spin-charge separation in the one-particle spectral 
function [68, 69]. The Bethe ansatz approach simplifies in the infinite U limit 
[70]. After earlier work [71, 72] the frequency dependent optical conductivity 
of the Id Hubbard model was also studied using Bethe ansatz methods [73, 74], 
as well as the dynamical density-matrix renormalization group [74]. 

5. Weakly coupled chains: the Luttinger to Fermi liquid 
transition 

Strictly one-dimensional systems are a theoretical idealization. Apart from 
this even the coupling to an experimental probe presents a nontrivial distur- 
bance of a Luttinger liquid. Unfortunately the weak coupling of a Id system 
to such a probe as well as the coupling between several LLs is theoretically 
not completely understood [26]. The coupling between the chains in a very 
anisotropic 3d compound generally, at low enough temperatures, leads to true 
long-range order. The order develops in the phase for which the algebraic 
decay of the correponding correlation function of the single chain LL is most 
slowly [64]. This can lead e.g. to charge-density wave (CDW), spin-density 
wave (SDW) order or superconductivity. 

In the following we shortly address some important issues of the coupled 
chain problem, which are a prerequisite for the theoretical descriptions of the 
attempts to experimentally verify LL behaviour. In the first part of this section 
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theoretical aspects of the problem of an infinite number of coupled chains are 
addressed. This is followed by a short discussion of the (approximate) experi- 
mental realizations of LLs. As there are other chapters in this book addressing 
this question the discussion will be rather short. 

5.1 Theoretical models 

We consider a systems of N± coupled chains described by the Hamiltonian 

N ± 

(ee) 

where the Hi are the Hamiltonians of the individual chains, the H-j repre- 
sent the two-body (Coulomb) interaction of electrons on different chains and 
the last term describes the hopping between the chains with t± i ij the 

transverse hopping matrix elements and the c|JL i the (creation) annihilation 
operators of one-particle states with quasi-momentum k n along the chain i 
and spin a (if spin is included in the model). The individual Hi can be TL- 
Hamiltonians Eq. (1.17) or lattice Hamiltonians like in Eqs. (1.63) or (1.66). 

We address the question if LL physics survives in such a model. The second 
and the third term on the rhs of Eq. (1.67) describe different types of cou- 
plings between the chains. If the transverse hopping is neglected (t± = 0) the 
model can be solved exactly for special assumptions about the two-body inter- 
action and the Hi. If the individual chains are described by TL-Hamiltonians 

(ee) 

Eq. (1.17) and the interaction iJ>- can be expressed in terms of the den- 
sities Pn t ( a ) t a,i the exact solution is possible by bosonization [75, 76]. This 
is important when the long range Coulomb interaction is taken into account. 
For a single chain the corresponding one-dimensional Fourier transform v(q) 
(which has to be regularized at short distances) has a logarithmic singularity 
for q — > 0. This leads to = and the divergence of the anomalous di- 
mension, i.e. the system is not a LL. The \kp harmonic of the density-density 
correlation function shows a very slow decay almost like in a Wigner crystal 
[77]. The Coulomb coupling between the chains removes this singularity and a 
three-dimensional extended system of coupled chains is a LL [75]. The corre- 
sponding anomalous dimension can be calculated and leads to values of order 
unity for realistic values of the coupling constant e 2 /(ttHvf) [76]. If 2fcp- 
scattering terms of the interaction are kept the model can no longer be solved 
exactly and a more complicated scenario emerges in the parquet approximation 
[78]. 

The inclusion of the transverse hopping presents a difficult problem even if 
the inter-chain two-body interactions are neglected. This is related to the fact 
that the transverse hopping is a relevant perturbation for < 1 [79-81]. This 
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can easily be seen if the individual chains are described by TL-Hamiltonions 
Eq. (1.17), scaling like 1/L. As discussed in section 3 the c~}\ . scale like 

(l/L) aL / 2 . As H( tj -} involves products of creation and annihilation operators 
on different chains no further boson normal ordering is necessary and 
scales as (l/L) aL . This suggests "confinement" for oll > 1: if an extra elec- 
tron is put on the j-th chain it stays there with probability close to 1 even in 
the long time limit. This conclusion can be questioned as RG calculations 
perturbative in t± demonstrate that the hopping term generates new and rele- 
vant interchain two-particle hoppings. These calculations show that the system 
flows to a strong-coupling fixed point which cannot be determined within the 
approach [81, 82]. 

If inter-chain two-body interactions are included the relevance of hopping 
terms can be different. When only density-density and current-current interac- 
tions between the wires are included, as discussed above [75, 76], the possible 
relevance around this Gaussian model, recently called sliding LL [83-85], can 
be different. If the single chains are in the spin-gapped Luther-Emery regime 
[57] single-particle hopping between the chains is irrelevant and the coupled 
system can show power-law correlations characteristic of a ld-LL [83, 85]. 
For the spinless model single particle and pair hoppings can be irrelevant for 
strong enough forward interactions [84]. 

In the following we concentrate on the Luttinger to Fermi liquid crossover. 
In order to get a quantitative picture it is desirable to study models which 
allow controlled approximations. The simple perturbative calculation in t± 
for the calculation of the one-particle Green's function by Wen [80] discussed 
below is unfortunately only controlled in the rather unphysical limit when the 
transverse hopping is independent of the distances of the chains (t± : ij = t±) 
[86]. The (retarded) one-particle Green's function G is expressed in terms of 
the selfenergy E 

g M-l.*;^) = - - i * - r — - y d- 68 ) 

where e, r denotes the energy dispersion for the noninteracting model and 

z = uj + iO is the frequency above the real axis. For small t± the dispersion 
can be linearized around k\\ = ±&f near the open noninteracting Fermi sur- 
face. This yields e fe ps ±vp{k\\ =F &f) + t±(k±). In the context of Fermi 
liquid theory the selfenergy is studied in (all orders) perturbation theory in the 
two-body interaction v around the noninteracting limit. This can be done us- 
ing standard Feynman diagrams. In the present context one wants to study how 
the LL behaviour for finite two body interaction and finite anomalous dimen- 
sion is modified by the transverse hopping. Similar to perturbation theory for 
the Hubbard model around the atomic limit nonstandard techniques have to be 
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used [87]. The simplest approximation, which corresponds to the "Hubbard 
I" approximation for the Hubbard model, is to replace £ in Eq. (1.68) in ze- 
roth order in t± by the selfenergy £( chain ) (/cy , z) of a single chain [80]. This 
approximation first used by Wen reads for kn ~ Uf 



G(k\\,k±, z; ijjwen 



1 



G+(fc||,. 



-t±(k A 



(1.69) 



where G + is determined by the spectral function p + discussed following Eq. 
(1.50) via a Hilbert transform. In the asymptotic low-energy regime this yields 
G + (kp+k\\, z) = A [(k^/k c ) 2 — (z /oj c ) 2 \ aL / 2 / (z— v c k\\) for spinless fermions, 
with uj c = k c v c and Aq = 7raL/[2sm (7tol/2)]. Wen's approximate Green's 
function leads to a spectral function with the same range of continua as p+(k\\ , u>) 
In addition there can be poles at ui k determined by setting the denomina- 
tor in Eq. (1.69) equal to zero. The poles located at u> k g ± = determine the 

Fermi surface k\\(k±) of the interacting coupled system. From Eq. (1.69) and 

the simple form of G + one obtains Ar J (&:||/A; c )( 1 ~ 0!L ) = t^(k±), which shows 
that the reduction of warping of the Fermi surface (FS) by the interaction is 
proportional to [^(/c^)/^] ^ 1-0 ^). This is shown in Fig. 1.6 for a two 
dimensional system of coupled chains. If one writes t±(k±) = t±c(k±), with 
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Figure 1.6. Fermi surface "flattening" in Wen's approximation for coupled chains for different 
values of the the anomalous dimension a l for a single chain. The dotted lines show the nonin- 
teracting FS, the long dashed curves correspond to ol = 0.125 and the full ones to az, = 0.6. 
At tti = 1 the FS degenerates to two parallel lines as without interchain coupling, called the 
"confinement transition". 



c(k±) a dimensionless function, the new effective low energy scale is given 
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by t c g = u c (t±/oj c ) 1 K 1 aL \ The weights Z^^ of the poles for k values on 

the Fermi surface are also proportinal to [t±(k±)/uj c ] aL ^ 1 ~ aL \ Wen's ap- 
proximate solution has the Fermi liquid type property of quasi-particle poles 
with nonzero weight on the Fermi surface, except at the special points where 
t±(k±) vanishes. The improved treatment by Arrigoni [88] shows that this 
peculiar vanishing of the quasi-particle weigths is an artefact of Wen's ap- 
proximation. The new idea involved is to let the number of "perpendicular" 
dimensions "d — 1" go to infinity. This extends the original idea of the "dy- 
namical mean field theory" (DMFT) [89], where one treats the Hubbard model 
in infinite dimensions as an effective impurity problem to the case of a chain 
embedded in an effective medium. Results are obtained by carrying out a re- 
summation of all diagrams in the t^-expansion which contribute in this large 
dimension limit [88]. This approach shows explicitly how the leading order 
Wen approximation is uncontrolled at low energies. For the case of weakly 
coupled one-dimensional Mott insulators one expects the approximation to be 
better controlled [90]. 

Despite the Fermi liquid like properties at energy scales much smaller than 
t e g the coupled chain system can nevertheless show LL like properties for en- 
ergy scales larger than t e g if there is a large enough energy window to the 
high energy cutoff uo c which describes the regime where the asymptotic LL 
power laws hold for a single chain. Then for temperatures lower than £u c but 
higher than t e s the system behaves like a LL. The integrated spectral functions 
p< (oj) probed by angular integrated photoemission, for example, show ap- 
proximate power law behaviour ~ (— Lj) aL for temperatures larger than t e g in 
the energy window ksT < —u < u c . Unfortunately little is known about the 
value of uj c for microscopic models. An exception is the Tomonaga model Eq. 
(1.14) with a constant v(k) up to the cutoff k c , where the high energy cutoff £o c 
equals uj c = mm(v c , v s )k c [4]. This implies for the integrated spectral func- 
tion for the very large U Hubbard model with periodic boundary conditions 
that the power law \u\ aL only holds in a narrow energy window ~ v s , which 
vanishes proportional to 1/U in the U — > oo limit [65]. Another example is 
the Hubbard model at boundaries where u c can be very small for small U [50]. 

As an alternative way to treat the "anisotropic large dimension model" [88] 
one can try to solve the resulting chain-DMFT equations numerically, using 
e.g. a quantum Monte Carlo algorithm [91]. In this reference the Hi were 
chosen as Hubbard Hamiltonians (1.66) with chain lengths of 16 and 32 sites. 
The results for a partly filled band as a function of temperature indicate in fact 
a crossover from a LL to a FL at the estimated crossover scale as the tem- 
perature is lowered. In agreement with Arrigoni [88] the authors find that the 
quasi-particle weight is more uniform along the Fermi surface than suggested 
by Wen's approximation Eq. (1.69). At half filling and low but finite temper- 
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atures the crossover from the Mott insulator to FL crossover was examined 
(the intermediate LL regime was too narrow to be visible). In the future it is 
to be expected that this method applied to longer chains and additional near- 
est neighbour interaction will provide important results which allow a more 
realistic comparison with experimental work. 

Because of space limitations the interesting field of a finite number of cou- 
pled chains cannot be discussed here [92]. 

5.2 On the experimental verification of LL behaviour 

There exist several types of experimental systems were a predominantly Id 
character can be hoped to lead to an (approximate) verification of the physics of 
Luttinger liquids. In the following we present a short list of the most promising 
systems and discuss some of the experimental techniques which have been 
used. As these topics are also discussed in other chapters of this book we do 
not attempt a complete list of references but only refer to most recent papers 
or to review articles on the subject. 
The following systems look promising: 

■ Highly anisotropic "quasi-one-dimensional" conductors 

There has been extensive work on organic conductors, like the Bechgaard 
salts [93, 94], as well as inorganic materials [95, 96]. 

■ Artificial quantum wires 

Two important types of realizations are quantum wires in semiconductor 
heterostructures [97, 98] or quantum wires on surface substrates [99, 100]. 

■ Carbon nanotubes 

The long cylindrical fullerenes called quantum nanotubes are also quantum 
wires but have been listed separately because of their special importance 
in future applications like "molecular electronics" [101, 102]. Using the 
peculiar band structure of the tt -electrons of a single graphite plane it was 
shown that single wall "armchair" nanotubes should show LL behaviour 
with K c ~ 0.2 — 0.3 down to very low temperatures [103, 104], despite the 
fact that two low energy channels are present. 

■ Fractional quantum Hall fluids 

Electrons at the edges of a two-dimensional fractional quantum Hall sys- 
tem can be described as a chiral Luttinger liquid [105]. The power law 
tunneling density of states observable in the tunneling current- voltage char- 
acteristics shows power laws of extraordinary quality [106]. The theoretical 
predictions for general filling factors between the Laughlin states v = 1 and 
v = 1/3 [107, 108] are not bourne out by experiment [109]. As in these 
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chiral LLs the right- and left-movers are spatially separated the edge state 
transport is quite different from the case of quantum wires and FQH fluids 
are not further discussed in the following. 

Promising experimental techniques to verify LL behaviour are: 

■ High resolution photoemission 

One of the earliest claims of possible verification of Luttinger liquid be- 
haviour was from angular integrated photoemission of the Bechgaard salt 
(TMTSF^PFg, which showed a power law supression at the chemical 
potential with an exponent of order 1 over an energy range of almost one 
eV [110]. There are serious doubts that this suppression can be simply ex- 
plained by the LL power law behaviour [94]. Therefore a large number of 
other quasi-one-dimensional conductors were examined [94-96, 111]. In 
addition periodic arrays of quantum wires on surface substrates were stud- 
ied by angular resolved photoemisssion (ARPES), but the interpretation of 
a two peak structure as spin-charge separation [99] was questioned [100]. 
Spin-charge separation was shown to occur in the Id Hubbard model also 
at higher energies on the scale of the conduction band width [69, 70, 73]. 
Recent ARPES spectra of TTF-TCNQ were interpreted with the Id Hub- 
bard model at finite doping to show signatures of spin-charge separation 
over an energy scale of the conduction band width. As for the Hubbard 
model K c > 1/2 for n / 1 which implies < 1/8 for the anomalous 
dimension the experimentally found nearly linear spectral onset at low en- 
ergies cannot be explained within the same model. ARPES data for the "Li 
purple bronze" seem to compare favorably to the LL lineshape [96]. For the 
quasi-one-dimensional antiferromagnetic insulators SrCu02 and S^CuOs 
ARPES spectra have been interpreted to show evidence of spin-charge sep- 
aration [112]. For a more in depth discussion see the chapter by Grioni in 
this book. 

■ Transport 

As discussed in section 3 even a single impurity has a drastic effect on 
the conductance of a LL, which vanishes as a power law with temperature. 
Another issue is the "conductance puzzle" of a clean LL. There has been 
an extended discussion whether the quantized value e 2 /h for noninteract- 
ing electrons in a single channel is modified by the interaction to K c e 2 /h 
[113, 114]. Apparently the answer depends sensitively on the assumptions 
made about the contacts, a very delicate theoretical as well as experimental 
problem [115]. Experimental results are available for cleaved edge over- 
growth quantum wires [97] as well as carbon nanotubes [116-118]. In the 
nanotubes the authors observe approximate power laws of the conductance 
which seem to be consistent with LL behaviour. A detailed dicussion of 
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transport through quantum wires is presented in the chapter by Yacoby. 
For a recent theoretical discussion of experimental results on the interchain 
transport in the Bechgaard salts see references [119, 120]. There the ques- 
tion of energy scales and the importance of the proximity of the incipient 
Mott insulator are addressed. 

■ Optical properties 

Optical properties have long been used to investigate electronic properties 
of quasi-one-dimensional systems [121]. The optical behaviour of differ- 
ent Bechgaard salts was analyzed recently using LL concepts [122]. At 
low energies, smaller than about ten times the Mott gap, the importance of 
dimerization and interchain hopping was pointed out [123]. As there is a 
separate chapter about the optical response in chains and ladders it will not 
be discussed further here. 

Obviously neither the list of systems nor that of methods is coming close to 
being complete. They were presented to show that there are intensive experi- 
mental activities in the attempt to verify the elegant LL concept put forward by 
theoreticians. Further work on both sides is necessary to come to unambiguous 
conclusions. 
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